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Abstract 

We study the Robinson- Trautman-Maxwell Fields in two closely 
related coordinate systems, the original Robinson- Trautman (RT) co- 
ordinates (in a more general context, often referred to as NU coordi- 
nates) and Bondi coordinates. In particular, we identify one of the 
RT variables as a velocity and then from the Bondi energy-momentum 
4- vector, we find kinematic expressions for the mass and momentum 
in terms of this velocity. From these kinematic expressions and the 
energy-momentum loss equation we obtain surprising equations of mo- 
tion for 'the center of mass' of the source where the motion takes place 
in the four-dimensional Poincare translation sub-group of the BMS 
group. 



1 



1 



Introduction 

One of the more important discoveries concerning integrating the Einstein 
equations was the realization that the condition of metrics being algebraically 
special allowed a huge simplification in the integration[l]. Not only was there 
the beautiful geometric result, the Goldberg-Sachs Theorem[2], saying that 
the degenerate principal null vectors must form a null geodesic congruence 
that is shear-free, but this result allowed all radial integrations to be per- 
formed. The algebraically special metrics separate into two classes, those 
with the null geodesic congruence having a vanishing[l] twist and those with 
non-vanishing[4, 5, 6] twist. We will study the subclass of the twist-free con- 
gruences, with a double degeneracy, known as the type II Robinson-Trautman 
(RT) metrics. Many of the issues explored here have been studied earlier by 
other authors[l, 3, 5, 6, 7, 8]. In particular Tafel and Drey looked at the 
relationship between the RT metrics and the Bondi metrics. 

It is the purpose of this note to describe some rather surprising results 
concerning the physical content that is contained in both the type II vacuum 
Robinson-Trautman[l] and Robinson- Trautman-Maxwell[5, 6] metrics, but 
had been essentially unobserved. Since, at least formally, the RT spaces are 
special cases of the RTM spaces, often we will just refer to both as RTM 
though there are apparently significant subtle issues that separate them. 

These metrics are largely determined by a single function of three vari- 
ables denoted hj V = V{t, C, C); where r is the NU/RT time[6, 9] coordinate 
and {(, () are complex stereographic coordinates on the sphere. (The coor- 
dinates used by Robinson-Trautman are a special case of those referred to 
as NU coordinates.) V{t,(,C)j whose geometric meaning will be clarified 
later, satisfies the Robinson-Trautman equation or its extension to include 
the Maxwell field. 

The basic plan in this work is to study the consequences of expanding V 
in spherical harmonics, i.e., assume that 

nT,c,c) = E^'"(^)^-(c,c), (1) 

and using the expansion in the RTM equations. The fact that the RTM 
equations are very non-linear means that the spherical harmonic decomposi- 
tion of the equation itself is quite complicated and that approximations were 
necessary. The approximated RTM equations for V then become a system of 
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coupled ODE's for the evolution of the components V'''^. The higher spherical 
harmonic components, i.e., those for / ^ 2, become derivatives of multiple 
moments. Our interest (partially) centers on the lowest components, i.e., 
(/ = 0,1), which can be interpreted as the four-velocity of a world-line in 
the four-dimensional space of the Poincare translation subgroup of the BMS 
group. The equations of motion for this world-line are contained in the lowest 
harmonic part of the RTM equations but are driven by the higher harmonics. 

The surprise comes when these equations are translated into Bondi co- 
ordinates and one looks at the Bondi energy-momentum four-vector and its 
evolution; namely the energy-momentum loss equation. Though the transfor- 
mation to the Bondi coordinates is complicated and can be done in practice 
only perturbatively in powers of v/c, nevertheless we see a beautiful confir- 
mation of the identification of the / = 1 components with the three- velocity 
V*. The Bondi three-momentum becomes, in the first order approximation 



In other words we obtain the usual kinematic definition of the momentum but 
now (by including the Maxwell field) with the standard radiation reaction 
term. In addition the Bondi energy, in the second order approximation, 
contains the relativistic kinematic mass-energy relation 



Higher moment contributions have been left out in both E and P\ 
The calculations that were involved in order to obtain these results were 
fairly lengthy and complicated and far from obvious. Since the equations 
are quite non-linear, considerable use had to be made of Clebsch-Gordon 
expansions of spherical harmonic products. In order to keep the presentation 
as clear as possible we have the following organization. In the second section 
we review our notation. In section three we derive the flat-space metric 
based on the null coordinates that are associated with the light-cones with 
apex on an arbitrary time-like world-line in Minkowski space. The variable 
^('T, C) C) appears naturally in this metric and one sees in this simple case 
how the velocity associated with the world-line enters the V{t, C, C) via the 
I — (0, 1) harmonics. In this flat space case the coefficients of the I = (0, 1) 
are essentially arbitrary while in our RT and RTM cases these coefficients 
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satisfy equations of motion and are driven by the higher harmonics. In section 
four we state the RTM metrics and their related differential equations and 
discuss their origin and how they are related to the general asymptotically 
flat Einstein-Maxwell metrics. In section five we state the details of our 
approximations and show how in this context the RT and RTM equations 
can be integrated. In section six we describe the transformation from the 
NU/RT coordinates to Bondi coordinates and construct the Bondi energy- 
momentum four vector and its evolution equation. Finally in the concluding 
section we discuss how this work fits into a much broader - e.g., to the twisting 
algebraically special type II metrics and even to a large class of algebraically 
general asymptotically flat Einstein-Maxwell fields. 

The calculations described here were often long and tedious and short- 
cuts were not found. This, unfortunately, does make the work difficult to 
follow in detail. To try to compensate for this we have attempted to give 
considerable background material. 

2 Notation 

Throughout this work we will be using certain Lorentzian and Euclidean 
vectors and tensor valued functions on the sphere. The most important of 
them is the null vector = r/"''/;, that sweeps out the entire sphere of null 
directions as 

(C,C)=co^^(e^^e-'^) (4) 
ranges over the sphere. Our space-time conventions are such that 

r]'''^r}ai,^diag{l, -1,-1,-1). 
In a given Lorentz frame we chose I"' to have the normalized form 

~ 2 1+cc' 1+cr " 

t"- and C" being unit time and space-like vectors with = (1,0,0,0). The 
factor ^ had been chosen historically so that the antipodal null vector n", 
with n"/a = 1, is 'symmetrically' defined, i.e., 

" - 2^^' i + cc' i + CC'i + CC^" 2^' 
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Remark 1. This normalization, which has been in use for many years 

and has many formal benefits, nevertheless has certain serious drawbacks, 
namely the appearance of some unpleasant factors of 2 and V2. Unfortu- 
nately they will be often seen and will create a serious annoyance. They 
appear, mainly, in the definitions of the several time coordinates. 

Note that for the covariant forms of l^ and n", we just change the signs 
in the last three components, e.g., 

,„ = ^(i,_i±4.iiz4,w£). (T) 

2 1 + CC' 1 + CC 1 + CC 

Other (co)vectors obtained directly from the la are 

3/ ^rn + -2C ^ 

"T.a = Ola — —^[0, =,— —, =), (8) 

2 ^ 1 + CC' 1 + CC 1 + CC ^ ^ 

^(1 + C^) -2C , 

'n^a = ola^—^iO, =, ^ =), (9) 

2 ^ 1 + CC' 1 + CC 1 + CC 

/2 

t- = ^(r + n«) = (l,0,0,0) = t„, (10) 
K, 1 + CC 1 + CC 

so that we can write 

, Ca , 1 1 , , 

We have that 



5r7(,) = Pt'ddPM^ ^Vis) = P^-^'d^{Po-\s)), Po = 1 + CC- 

Note that in the given Lorentzian frame, we can treat nia, fria, and Ca as 
Euclidean three- vectors which we write, with Euclidean indices i,j , as 

rui, Wii, and q. 

Also note that from the definitions of the stereographic coordinates, Eq.(4), 
we have that 

Ci — k — rii — — -\/2(cos (j) sin 9, sin sin 6, cos 6) — —V2C^ 
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where C" is the unit Euclidean radial vector. 

The unifying notation of the tensor spin-s harmonics [10], Yi\j sYim 

which allows us to write 

Yo' = 1, (13) 
Yii = m, (14) 
y° = QY,]^dY,-'^-c,, (15) 
Yu' = nii, (16) 

and 





— mirrij, 




(17) 






-{cirrij + miCj), 


(18) 






3CiCj Qidijf 


(19) 




= -{cifrij 


+ miCj), 


(20) 


y-2 


— mifrij, 




(21) 



turns out to be almost indispensable. The ordinary spherical harmonics Yi^ 
are equivalent, i.e., are linear combinations of the s = 0, Yilj j^, 

Ylm <^ ^ij..fe 

SO that, for example 

Ylm = (1^1^10,^1-1) ^ Yii = -a- (22) 

Due to the non-linearity of the RTM equations, products of these har- 
monics often appear and must be reduced via a Clebsch-Gordon expansion. 
Several of the most important and most used of these expansions[10] are: 



yly-l 



—dj. 
3 - 

\ + -YL 
3 3 



2ii' 



(23) 
(24) 



and 
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^2W^27/ = -^{2^ik^ji + 2^ii^kj - i^SijSki} (25) 

+ {I — 3 and 4 harmonics). 

Using the tensor spin-s harmonics, the spin-s functions can be expanded 
via the y^^^ ^. In particular, for the ordinary functions on the sphere, i.e., for 
s = functions, we have that 

F(o) = /°ro° + + r >^2°,- + 

rather than the equivalent expansion in terms of the ordinary spherical har- 
monics. A particular function that often appears is 

G(o) = gX = -^g' - \g'yl (26) 

where we have used Eqs.(12) and (15). 

We will be using several different time variables (or coordinates) in the 
course of this work, i.e., r, ub- The variable r is the proper time along 
a four-dimensional world-line; its derivative is denoted by 9^ = ('). Often we 
will be dealing with functions of w and of r that functionally are identical, 
i.e., F{t) and F{w) are identical functions. In that case we will again use 
the same notation, 9^ = (') with the meaning coming from the context. The 
(rescaled) Bondi time, ub — w/\/2 will be denoted by d^g — ( )■ The origin 
of the-\/2 is related to remark 1. 

3 A Flat Metric 

Before considering the RTM metrics we describe a flat-space version that 
imitates many of the features of the RTM metrics, albeit with a great simpli- 
fication. Starting with Minkowski space and ordinary Minkowski coordinates, 
x", we introduce a coordinate transformation to new coordinates (r, r, C, C) 
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= + + (28) 



that are based on the null cones emanating from an arbitrary time-like world 
line which we •parametrize so that it has a unit velocity vector, 

V^Va = 1, (27) 

so that (with our approximation, i.e., for slow motion) 

1 . . Iv'^ 

^O(r) = J v^'dr = T + S^". (29) 
The transformation is given by 

x« = r(T) + fr(c,c), (30) 

with x"' — ^"(t), the parametrized world-line, /"(CC); the null vector of 
Eq.(5) and Vo{tX,C) defined as the r dependent spin- weight s = function 
on the sphere 

Vo = C'la = vX, (31) 

r' = ^ 

^ dr ' 

Remark 2. If the world-line was simply the time axis, i.e., ^^{ub) — 
UBt"', then the transformation, x"" <^ (rts,r, C,C) ^■^ given by 

= MBt"+ V2fr(C,?). (32) 

In this special case, the new coordinates are the (flat-space) Bondi coordi- 
nates. 

To construct the new form of the metric we first take the differentials of 
Eq.(30) and simplify it: 

dx^ = i^'dr + d{^)nU) + {rdnco- 

Using 

c^r = PQ\m''dC + rff'dC), 
dVo = V^dr + P^^v''{madC + madC), 
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from Eqs.(8) and (9), we have that 



dx"" = v'^dr + l'^i—- -^rdr) (33) 



Tv\mkdC + mbdC) + ——{rrfdC + nfdC), 

which is then substituted into ds^ = rjabdx'^dx'^ . Using the scalar products 
v°/a = Vo, m"'rna = 1 and m"/a = we obtain [after the cancelation of two 
terms both containing v'^{mbd( + TfibdQ] the flat metric 

ds"" = r]abdx''dx^ = (1 - 2^r)dT^ + 2dTdr - r'?^, (34) 

or (with P = ^) 

ds" = r]abdx''dx^ = (1 - 2^r)dT^ + 2(iTdr - r'^. (35) 

The Gaussian curvature of the r = constant surfaces is, after a brief calcu- 
lation, 

K = iP^d^d^logP = SaiogPo = 1, 

the curvature of the unit sphere metric. 

Remark 3. It should be emphasized that in the present flat-space case 
we have the expressions: 

Xo = C{T)la = ^{r + 5e) - ICYl (36) 

Vo = v'^{T% = ^il + Sv')-lv%l (37) 

= v^'{r)k^^Sv"-lv^'Yl (38) 

They will be generalized to include higher harmonics in the next section. 
The metric form will appear quite similar. 

To make this metric, (35), agree with a certain conventional form[6] 
(coming from the normalization) that has been frequently used, we make 
the following simple rescaling of the r and r 
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resulting in 



with 



r = ^&r = \/2r*, (39) 
v2 

a _ 

V - ^° 



V*' r* 2 

ds^ = 2(/r--^r*)(ir*2 + 2(ir*dr*-— , (40) 

X = AP^ d^d^log P ^ I, (42) 



Note that we continue to use (') for both r and r* derivatives. The 
distinction should be clear from the context. 

As we remarked earlier for the flat case, our function Vq = v"(r)/a will be 
generalized in the RTM case to include higher harmonics. In the flat case, V 
codes the information of the velocity vector of the Minkowski space world- 
line. In the RTM case there is no such 'real world-line' but nevertheless a 
'world-line' in an auxiliary space does exist. It is easiest to show its meaning 
first in the flat space case. To do so we first return to the fiat space Bondi 

coordinates (MB,r,C,C),. Eq.(32). Rescaling r and ub as in Eq.(39) 

x" = v^u*Bt'' + f*/"(C,f). (43) 
The relationship between {u*Q,r*,(X) {'^■>i^->CX) is given by 

which with a bit of effort can actually be disentangled as a coordinate trans- 
formation. We need it only at future null infinity, J"*", i.e., in the limit as 

r^oo, where ((,() ^ (CO- 
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In this limit, if (44) is multiplied by la we obtain the relationship 

u*B = e{r)la{CO. (45) 

But since the 4-d Poincare translation subgroup of the BMS group is given 
by 

we see immediately that we can interpret the 'world-line' as a curve in the 
space of BMS translations, d". 

We will see how this type of 'world-line' arises naturally in the study of 
the RTM spaces. 

Remark 4. Often there are two different time- coordinates that are used 
at null infinity, 3^, namely, the Bondi time Ub and the so-called NU time r. 
Roughly speaking, the 2-surfaces of constant ub are uniformly stacked above 
each other while the constant r surfaces have an arbitrary uneven stacking. 
They are related by an expression of the form, ub = X{t,(,0- Eq.(45) is an 
example. This description applies to both fiat and asymptotically fiat space- 
times. 



4 RT and RTMaxwell Space-Times 

4.1 The metric and differential equations 

The RTMaxwell metric, after all the radial integrations have been performed, 
can be put into the form[l, 5, 6], (see Eq.(40)) 



ds'' 
with 



V*' , V*^'ibf 



+ 



v*^(t)i Vi 



)dT*' + 2dT*dr*- 



r*' dCdC 



, (46) 



K = AP^d^d^logP, (47) 

P - (48) 
= -V*-'x, (49) 



and Maxwell field 
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0S = 0, 

(f)l = 0f = |^*-V-^ (50) 

The two gravitational variables, V* and x = — T^*^'V'2°' the Maxwell 
field, 02°! functions only of (t*,C,C): ^^id satisfy the three coupled dif- 
ferential equations, the RTM equations, 

0t°' + 5(.*)[F>;°] = 0(51) 
+ 2fc\^*3^f = 0(52) 

{v*-'xy - - y*-%*)y* ■ ^r*)V*} + = 0(53) 

The constant q, in Eq.(50), is the conserved electric charge while x (or V'2°) 
carries the information of the Bondi energy-momentum four-vector which 
will be extracted later, k = 2G/c'^, and 5(^.) is the edth operator holding r* 
constant. The quantities and asymptotic tetrad components of 

the Maxwell field while x (or '02°) ^^e asymptotic tetrad components of the 
Weyl tensor. 

It is useful to note the logic of these equations: They can be solved 
(in principle) individually, in order. The first allows (fif to be expressed 
as a function of V while the second does the same for x- Finally, the last 
determines the r-evolution of V. 

Remark 5. The RTM equations are extremely difficult equations to solve 
or analyze rigorously. Very few exact solutions have been found. In the special 
case when q = 0, i.e., the pure RT equations, Piotr Chruscielfll] showed 
that, with sufficiently 'nice' initial data, the solutions evolve asymptotically, 
in time, to the Schwarzschild metric. In the general case, though one would 
hope for a theorem saying the evolution would lead to the Reissner- Nordstrom 
solution, this apparently has not yet been shown. It does make an attractive 
conjecture. Later, in our discussion of the approximate solutions, we can see 
how this issue arises. 

Remark 6. For the time being, to streamline the notation and avoid a 
plethora of similar symbols, we will drop the (*) on the V and r. Later, to 
see the physical relevance of our results, we must remember to restore the 
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original V and r via Eq.(39). In addition, the velocity of light, 'c', will be 
explicitly introduced via t ^ ct and w =^ cw. 



4.2 Background and Derivation 

In the remainder of this section we will attempt to place the RTM equations in 
a broader context, i.e., study how they are related to general asymptotically 
flat Einstein-Maxwell fields and how from that point of view they can be de- 
rived. We begin with the general asymptotic version of the Einstein-Maxwell 
fields in a Bondi coordinate and tetrad system and see how to specialize to 
the RTM equations. Later, using these results, we will reverse the procedure 
and go from the RTM equations to the Bondi point of view in order to obtain 
the physical meaning of the RTM variables. 

4.2.1 Bondi coordinates and tetrad 

We start with the full set of asymptotic Maxwell and Bianchi identities in 
Bondi coordinates, {ubX)C) and Bondi tetrad n, m, m)[6, 9], (the ub is 
already rescaled) 

Maxwell Equations 



- cj^l = 0, 
= 0. 



(54) 
(55) 



Bianchi Identities 



da', 



—a 



(56) 
(57) 
(58) 
(59) 
(60) 
(61) 



^2° - 



S Cr — 5^0" -|- (7 (7 — (7 (7, 



where 




(62) 
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From the definition of the mass aspect 

^ = ^0 +g2^+aa, (63) 
we see that (61) is equivalent to 

^ = (64) 

and by eliminating V''2 iii terms of ^, Eq.(58) becomes the very important 
Bondi mass loss equation 

*- = (j-(a)- + A;0°^°. (65) 
4.2.2 Asymptotically shear-free tetrad 

These relations take a new form after using the tetrad transformation, (the 
asymptotic null rotation around n), 

r = r + 6m" + 6m" + (66) 

h = _L/r + 0(r-2), 
L = L{ub,C,0, 

with 



0° = ^l' + 2L<j)l' + L'cl>l', (67) 

0? = + (68) 

00 = 0;°, (69) 

= Vo° + 4i^^f + 6i^'^2° + 4i^'^3° + ^Vf , (70) 

^0 = + SLT/^f + 3LV3° + ^^V'f , (71) 

= + 2LV;;o + LVf , (72) 

V'" = V^3° + ^<, (73) 

V'4° = <■ (74) 



Remark 7. A^'o^e that the ( *) used here has nothing to do with the ( *) 
from the rescaling of V and r in section 3. 
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First, we obtain a new asymptotic shear[12, 13] a* for the null vector T", 
that is given by 

a* = cj- [dL + LL). 

The requirement of a vanishing a* leads to the differential equation for the 
determination of L, 

SL + LL = a. (75) 

Using such an L, with Eq.(66), the asymptotic Maxwell and Bianchi Iden- 
tities are changed to equivalent equations (just with a changed tetrad) but 
with the same Bondi coordinates. Now however defines an asymptoti- 
cally shear-free and, in general twisting, null vector. Shortly we consider a 
special case where the twist vanishes. 
Maxwell Equations 

(0*°)- + 9(/.f + 2L-(/.f + L«)- = 0, (76) 
(0f)- + g0;o + (L0;°)- = 0. (77) 

Bianchi Identities 





= -5<-L«)- 






--*o 


(78) 


«)■ 


= -g^f -L(^;0)- 


- 3i:-^2*o 


+ 2k(j)f< 


--*0 

92 , 


(79) 










-*0 
'2 , 


(80) 




= d{ay+L{ay, 








(81) 












(82) 




= (5 a — S^a + aa 


— aa 






(83) 



+2Ld{ay - 2L 5(a)- + lV" - LV', 



where now we treat L{ub, C) the basic asymptotic data and a is defined 

by 

a = QL + LL. (84) 
By using (72) in (63) we find that 

^ = + 2L5(a)' + LV' + + aa, (85) 

so that (83) is again equivalent to 

^-^ = 0. (86) 
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Note also that only the first three of Eqs.(78)-(82) have a dynamic content, 
Eqs.(81) and (82) are just definitions that can be inserted into Eq.(80). This 
fact will be used later. 

4.2.3 Type II metrics in Bondi coordinates 

Up to now the two sets of equations, with the Bondi tetrad and the shear-free 
tetrad, are completely equivalent. We now make our major condition that 
is a severe restriction on the solutions. We assume that we are dealing with 
algebraically special type II Einstein- Maxwell metrics with a non-twisting 



pnv, This translates immediately into the statement that V-'o" = V-'i^ = 
and that automatically, from Eq. (78) , we have that 0^° = 0. The non-twisting 
assumption has not yet been used. 

Using this the Maxwell/Bianchi Identities become 

g(/)f + 2L-0f + L(0f)- = 0, (87) 
«)- + S0;° + (L0;O)- = 0, (88) 

ir^'y = -Qi^;'-L{i;;'y -2Lrs'+k<i>i%\ _ (90) 

* = V'2° + 2Lg(a)- + L'(a)" + g'a + (7(a)- = (91) 
with, as mentioned above, the defining relations 

a = dL + LL. (93) 



They are a set of coupled equations for the variables (L, 0^°, 02°; '02°)- 
The shear is no longer free data - it follows from knowledge of L, via 
Eq.(93). 

4.2.4 Transformation to NU/RT coordinates 

By transforming to a new set of coordinates at 3'^, namely NU/RT coor- 
dinates, [with the same tetrad], we obtain a great simplification in these 
equations. 
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The transformation is to be given by[13, 9] 

ub = X{tXX), (94) 

{( unchanged) where X(r, C,C) is obtained from the following construction, 
based on the (assumed) known L{ub, C^O '■ First the (CR[14]) function, 

r^T{uBXX): (95a) 

is found by solving the differential equation 

QT + LT=0, (96) 

and then by inverting it, i.e., by algebraically solving for Ub, the X(rX,C) 
is constructed. 

Remark 8. In the more general situation of twisting type II congruences, 
the T and X are in general complex-valued analytic functions, however for 
the (non-twisting) Robinson- Trautman case they are real. 

Now applying this transformation, ub =^ t, to Eqs.(87)-(91), (a lengthy 
process), we obtain our final version of the Maxwell/Bianchi Identities in 
the new tetrad basis with the NU/RT coordinates, (r, C,^). The variable 
L{ub,CX) is replaced by the new basic variable 

V{T,C,O^X'{r,C,0- (97) 

Remark 9. Frequently, in the application of this transformation, we have 
used the easily derived relations 

dT 

L = -_ = g(,)X, (98) 

L = d^r)V, (99) 

a = dL + LL = 5^^)X, (100) 

where d(j^ means the d- operator hut holding r constant as opposed to holding 
Ub constant. 

4.2.5 Maxwell/RT in NU/RT coordinates: 

^ir)[V^f] = O^0f (101) 
cPf + Q^r)[VcP;'] = 0, (102) 
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-S(r)X 

X 



2A;\/30f 02 , (103) 
(104) 

(105) 



with (x, V, t/'2°) i^sal. The mass aspect simplifies, after a lengthy calcula- 
tion, to 



These equations are precisely the same as those given earlier, Eqs.(51 
)-(53). They had been obtained[6], with some notation changes, by direct 
integration of the NP equations assuming that V^q^ = — — with a 
real 'divergence' p. 

4.3 Approximate solutions 

As we mentioned (and is easily seen) the RTM equations are quite non- 
linear and difficult to solve. We have resorted to an approximation scheme 
based on the idea that the dominant term in gravitational radiation is from 
the quadrupole. We thus assume that all quantities are to be expanded in 
spherical harmonics including terms from / = 0,1,2 with nothing higher. 
Essentially we are approximating ofi" the Reissner-Xordstom metric which we 
take as zeroth order. We consider terms up to quadratic in the deviation 
from the zeroth order. Due to the non-linearities there are frequent products 
of the I — 0,1,2 harmonics which are then decomposed via the Clebsch- 
Gordon expansions (Sec. II). When these expansions involve harmonics with 
/ > 2 these higher / terms will be omitted. In a few instances where the 
/ = 2 harmonics or the quadratic terms do not appear to be important, we 
will omit them. This occurs for the Maxwell fields, which appear in the 
gravitational equations only as quadratic terms. These I — 2 could easily 
be restored but to keep the equations tractable their omission seems to be 
reasonable. Where this occurs, we will point it out. 

We use the following notation (remembering the rescaling of r) for the 



(106) 
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harmonic expansions of the relevant variables: 



^ = ^ - \vX^ + v-n, = ^^^^ - + .-y4, (107) 

V = ^-^-Iv^n + v'^'Yir (108) 



(109) 



e = ^^ = (no) 

v<f>*' = ^^ji + ^^^y^^S (111) 

X = Xo + <^Xo + x'l^i° + x''>^2°r (112) 
Our first difi'erential equation, Eq.(51), 

<Pf + Q^r)[V<Pl']^0, (113) 

becomes, using the properties of edth, 

-ilv-'y = g^-v = a(,)[i-0;°] = S(r)[*Y-i + *^^i^2^], (114) 

From the Taylor expansion of V~^, with (108), we have the solutions, 
V<f>l' = -q[V2v"-^{v^v'^y]Y-' (115) 



+q !^2V2v^'' - Sym.T.F.[vh'' + i^lKl^^^'^^^^'j | y-/. 

Sym.T.F. means symmetrized and trace free part of the expression and xo 
is constant. 

The linearized Maxwell field is 



C = 0, (116) 
6f = g[l + v^t;%°-2v^t;^%°J, (117) 
bf = -2q[v%.' - 2v^^%.}]. (118) 
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(See appendix A for the solution up to second order with the I — 2 
harmonics.) 

Using Eqs.(50), (112) and (115) in Eq.(52), i.e., 



g(,)X + 2A;y V02° = 0, (119) 



we find that 



X' = -^(x/2^^'-y{^Vn'), 

X'^ = '^f^v^^'-t^Sym.T.F.{v^v'' + —vV). 
3 6 7 

Finally the last and most complicated equation, (53), is studied after 
inserting ^^id x from Eqs.(115) and (112). We find, after considerable 
effort, that the evolution is given by 



(<^Xo)' = SxovV' + ^^v'^v'^ (120) 

o o o 

^^yi' = -^^yjyij +q'^k{-[v"' -^^{V^v'^y (121) 

5 15 

48 • 30 kn'^ 
3Xov''' = -6v'^ - Sym.T.F.[^^v%^^] + ^v''" (122) 

7v 2 3 

+Sym.T.F.[^^^(vV - 3vV - —v%^' - —v^%'% 
^ ^ 7^2 64 2 

Equations (120), (121) and (122) constitute our final set of (ordinary) dif- 
ferential equations within our approximation scheme for the RTM equations, 
while Eqs.(112) with (116) and (50) yield the asymptotic physical Maxwell 
and Weyl fields. Even with our 'heavy' approximations it is not easy to un- 
derstand them. However, in the next section, where we transform the RTM 
equations to Bondi coordinates, we can see meaning to these equations from 
physical insight. This physical insight might even help in understanding the 
mathematical structures involved. 
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5 The Bondi Energy-Momentum Four- Vector 



In order to construct the Bondi energy-momentum four-vector and its evo- 
lution we reverse the coordinate transformation that took us from Bondi 
coordinates to the NU/RT coordinates that was given earlier in Eqs.(94) and 
(95a). We then express the mass aspect in terms of the RTM variables. 

5.1 Transformation to Bondi coordinates 

In order to go from the present NU/RT coordinates, (r, (, Q on null infinity, 
(3"*"), to the Bondi coordinates, {ub,CX)j simply use Eq.(94) with its 
approximation 



UB = X{t,C,0 (123) 
= f V{r'X,Odr' ^^-^-e{r)Y,', + e'{r)Yl, (124) 

with (CO unchanged. We need to find the inversion of (124), i.e., we must 
find the approximate r = T{ub, C) 0- 

Remark 10. // exact solutions, V{rX,C): of the RTM equations were 
known, then the exact transformation to Bondi coordinates at 3^ would he 
given by 

ub = x{tXX) = fv{T'XX)dT', 

with an additive arbitrary supertranslation. By ignoring it, we implicitly have 
chosen a specific Bondi frame. 

Remark 11. In this section we will be returning to the unrescaled r, so 
that T ^ t/V2, (see Eq.(39), (107) and (109)). 

Before finding T{ub-iC-iC)-i point out that OubT = T' — since by 
differentiating (123) 

l^V-^^VT-. (125) 

dUB 

The relation r = T{ub, C,C) can be approximated in the following man- 
ner: from Eq.(107), we have that 

e{r)^r + 5e{r), (126) 
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F(r,C,C) ^ -5e{r) + ^C{r)Y^i-^/^CKr)Y^ir 



T = nj-5e{w) + ^C{w)Y?i-^CK^)Y^ip (131) 



so Eq.(123) can be rewritten, with 

w = V2ub, (127) 

as 

T = «; + F(r,C,C), (128) 
2 

This can be iterated with 

T^w (129) 
as the zeroth order iterate. The first iterate of Eq.(128) is then 

T = w + F(w,C,C), (130) 

V2. 
2 

with the second iterate 

T = w + F{w + F{wXX)XX) ~ w + F{wXX) + Fd^F, (132) 

with higher iterates (though unpleasant and not needed) easily found. We 
then have that 

T ■ = ^d^T = V2(l + d^F{w) + dlF{w) ■ F{w) + d^F{w) ■ d^F{w)), 
or, using the Clebsch-Gordon expansions and omitting the argument w, 

T- = V2{1 + -{v'CY -Svo + -^[C'v'^' + v'^v'^] (133) 
3 5 

2 5 

6 7 

This result plays an important role in determining the mass aspect and 
the Bondi four-vector. 

Remark 12. We have been using the notation (') = dr and (') = duB- In 
the following where w, from Eq.(127), is used in the same functional expres- 
sions that had used r, we will continue to use the notation 
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5.2 The Mass Aspect 

Our task is to take the mass aspect from Eq.(106) 

qr = ^ = r2' + ^r)%)X (134) 

which is a function of (r, (, Q and reexpress it, via r = T{{ub, C, C); in Bondi 
coordinates. 

All the terms on the right side of (134) are known functions of (xo) °) v'', v^^ , 5^°, 
^*-^), which in turn are functions of (r, These known functions are sub- 
stituted into Eq.(134) and then r is replaced by the first iterate, Eq.(131). 
{The first iterate is sufficient since all the relevant quantities are already 
order.} 

After another very lengthy calculation, using the third power of T ', i.e., 
Eq.(133), and frequent use of Clebsch-Gordon expansions, we finally find 
(including only the zeroth and first harmonic terms, the only ones needed for 
the Bondi four-vector) 

*(w,C,C) = ¥'^ + ¥'^^Y,l (135) 

1 144 288 

+5X0 - q'k[^v'"e + fc'v'"' + fv-'v"^' + v'v']}, 



3'v^ 5' ' 5^ 5 ^ 5 

Remark 13. As already pointed out the already defined w = \/2ub is 
really the unrescaled time coordinate, Eq.(43), that was originally used. In 
what follows we replace w by cw, i.e., 



w =^ cw, 

C) ^ c-H') 



Since the relationship [9, 15] between the mass aspect and the Bondi 
energy-momentum four-vector, (Mc, P') is given (with M — Mq + SMq and 
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Mo = const.) by 



^(w, C, C) = - ^P%^ + (136) 



then by comparison with Eq.(135), we can make the physical identifications: 
with 



G 



we have that 



M = Mo[l + — + ^ + — (eV^-' + 2^^V^)] + ^eV^- (137) 



,2 



3(f be 5G 



+^-4-[12(v^t;^^)' + 3f + + 6v*%^']. 



5.3 The Bondi Energy-Momentum 4- Vector 

The Bondi energy-momentum four-vector, 

P» = (Mc,PO> (139) 

which has been extracted from the asymptotic Weyl tensor has been known 
for almost 50 years. It transforms as a Lorentzian four-vector under the 
Lorentz subgroup of the BMS group. It satisfies the beautiful Bondi energy- 
momentum loss theorem. But, to our knowledge, it never before had a kine- 
matic relationship with anything that resembled a velocity. 

We see that the variable that was introduced formally, in Sec. IV, as 
the I — 1 harmonic component in V{t, Q now plays the kinematic role of 
a velocity. We see that the mass has the kinetic energy term 



Moc\l + --)^^^, (140) 
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while the momentum, at lowest order, is just 



P' = Mqv\ (141) 

Though most of the remaining terms in both M and come from the 
I = 2 (quadrupole) terms and remain to be understood, the several terms 
involving the / = 1 quantities (both linear and quadratic) are interesting. In 
M there is the anomalous term 

which can be interpreted as Force x displacement, the work done on the 
total system moving it to the 'position' ^. The two quadratic terms appear 
to be new, a contribution to the mass from an accelerated massive charge. 
This appears to be a genuine prediction from this work, though the effect is 
probably too small ever to be measured. 
The second term in P\ i.e., 

-^?V', (142) 

is quite well known. In the equations of motion, obtained below, from the 
energy-momentum loss theorem, the derivative of this term is precisely the 
(electromagnetic) radiation reaction force with the correct numerical factors. 

The question is in what 'space' is this motion taking place? We saw in 
Eq.(124), i.e.. 



UB = X{r,C,0^ (143) 



V2 



that we omitted an arbitrary additive (integration) term, coming from the 
four-parameter Poincare subgroup of the BMS group, of the form 

Cola = ^ - (144) 

where is an arbitrary (constant) translation. We can thus interpret the 
^"(r) as motion in this four dimensional space and define it as the center of 



25 



mass motion of the total system with v"(t) = C"'(t) as the center of mass 
velocity. 

Remark 14. Though this issue will not be pursued here, there is an 
alternate space where the motion or world-line can be considered as taking 
place. Associated with each RTMaxwell solution there is defined a complex 
four-dimensional manifold referred to as H-space[13, 16, 17]. Our world-line 
can be interpreted as a curve in H-space. 



5.4 Energy-Momentum Loss Theorem 

In a Bondi coordinate and tetrad system one can find the Bondi energy- 
momentum loss equations directly from the asymptotic Weyl tensor equation[15, 
9], Eq.(65), 

= a-{ay + k4'^^, (145) 
k = 2G/c\ 

By substituting, from Eq.(118), 
and, from Eq.(lOO), 



a = V224—Y' (146) 



c 



then using, on the left side, Eq.(136), 



^{w, C, C) = - ^ + (147) 

we have, (after Clebsch-Gordon expansions, from Eqs.(23) and (25)), the 
mass (energy) loss theorem 

M' = 5M' = -^cv'^v'^ - \q^v"v" - ^v'^'v'^', (148) 
5G ?>& b& 

or 

E' = {Mc^y = -'^^^v'H'^ - %v''v" - ^v'^'v'^', (149) 
5G Sc"^ 5c'' 

and the momentum loss equation 
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pi> ^ ?^L^yi'yi3\ (150) 



8\/2g' 
5c' 

The c has been explicitly introduced. 

There are several important things to notice about these relations: 



• the energy loss due to the electromagnetic radiation is precisely the 
electric dipole radiation expression, [with the dipole moment being 

d' = qC, (151) 

allowing us to identify as the center of charge as well as center of mass] 
and electric quadrupole radiation, if we identify the electric quadrupole 
moment, D^^ , with our from 

so that[18] 

— '—v'-^v^^ = = — 

5c3 180c5 

the classical quadrupole electromagnetic energy loss. 

• By comparing our term, from Eq.(149), for the gravitational energy 
loss, namely 

288 Q • • • 

c v^v^, 

5G 

with the standard quadrupole radiation formula[18], i.e., 

^,{Q''"W"}, (152) 
oc 

we find the physical meaning to the variables, namely modulo nu- 
merical factors, it is the 2'*'^ time derivative of the mass quadrupole 
moment: 

^ Q'^", (153) 



,4 



12\/2c' 
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Remark 15. Note that both the gravitational quadrupole moment Q^^ 
and the electric quadrupole moment D^^ are defined from the RT variable 

The equations of motion for the world-line arise when P*, from Eq.(138), 
is used in the momentum loss equation. Writing Eq.(138) as 

P' ^ Mv' - ^qh'' + MG + q^H, (154) 

where Mq has been replaced by M and the G and H contain I — 2 terms 
(which will be neglected in this discussion), P*' = ^^^-v^'v^^' becomes 

Mv'' = \q\'" + ^^:^v^'v'i' - M'v' (155) 

Using the mass loss equation, Eq.(148), it becomes the equations of motion 
for the world-line 

3c"^ 5c* 3c^ 5^^ 5G 

(156) 

Many of the terms in this equation are quite familiar from classical elec- 
trodynamics. If we ignore the mass loss terms (the three cubic terms) we 
have the well-known radiation reaction equations[18] whose solutions have 
the troubling run-away exponential behavior. By adding in the just electro- 
magnetic radiation term (the first and second cubic term) we again have a 
classical electrodynamic result [19] which was derived in the past by model 
building and which required mass renormalization. This term has the cor- 
rect sign and structure to possibly suppress the rim-away behavior - though 
whether or not it does so is not clear. The gravitational radiation term (the 
third cubic term), arising from electromagnetic quadrupole radiation, simply 
add to the likelihood that the run-away behavior could be suppressed by 
these equations. At this point it is not clear what effect the quadratic term, 
which couples the / = 1 and / = 2 harmonics, has on the motion. 

There are two important comments to make concerning Eq.(156), one 
positive, the other negative: 

• The negative comment is that we have cheated a bit, in that we have 
included a cubic term in Eq.(156) that comes from the mass loss, while 
throughout this work we have excluded cubic terms. Our defense is 
that it arises naturally from the mass loss and has such a clear physical 
meaning. 
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• The positive remark is that we have not used any model building, as 

has been used in the classical electrodynamic derivations, with the often 
used mass renormalization. We have here simply taken the Einstein- 
Maxwell equations and applied them to the class of non-twisting alge- 
braically special metrics and obtained the classical equations with the 
additional gravitational correction term. All the numerical factors are 
identical to those of the model dependent derivations. We have nothing 
that even resembles mass renormalization. Our mass is the Bondi mass 
seen at infinity. 

It has been known that there have been difficulties[ll], as we mentioned 
earlier, in giving a definite answer to the question of the end state of an 
RTMaxwell solution with a non-vanishing charge. It appears likely that 
this difficulty is associated with the delicate balance between the run-away 
behavior and its suppression via the three radiation terms. We hope to be 
able to study the numerical evolution of Eq.(lo6) in the near future. From a 
physical point of view we know that charged particles do not self-accelerate 
so there is some reason to believe that the same would be true for this special 
class of fields. 

6 Conclusion 

In this work we have returned to the study of the beautiful class of metrics 
known as the Robinson-Trautman and Robinson- Trautman-Maxwell equa- 
tions and tried to give them physical content. Many attempts over the past 
forty or so years have been made to solve them - with very few exact solu- 
tions - most work being on approximations. Existence theorems for the pure 
vacuum case were given and much approximate work was done on them but 
there appears to be little progress on the Maxwell version. There have been 
several attempts[7, 8], via approximations, to relate the pure RT metrics to 
the more general class of Bondi metrics where one has an understanding of 
asymptotic energy and momentum as well as the multipole structure. How- 
ever, to our knowledge, little or nothing along these lines was done for the 
RTM equations. 

In this work we have first studied the RTM equations, in their natural 
NU/RT coordinates, and reduced them, using a truncated spherical har- 
monic expansion and approximation, to a series of ODE's. Afterwards we 



29 



applied a method, recently developed, to approximate the transformation to 
Bondi coordinates and finally obtain the Bondi energy-momentum 4-vector 
and its evolution. By studying these equations, we were able to give physical 
meaning to many of the RTM variables. In addition we obtained directly 
from the field equations, with no model building, the classical particle equa- 
tions of motion with the radiation reaction force and terms for the possible 
suppression of the run-away behavior. 

There are several different directions in which we intend to extend this 
work: 

• We would like to redo the calculations done here consistently to third 
order mv/c. The obstacle is the complexity of the calculations. In the 
present work often we had to do the calculation two or three times to 
get agreement. 

• We plan to generalize this work to type II algebraically special met- 
rics that now contain shear-free null geodesic congruences with non- 
vanishing twist. It appears virtually certain from earlier work and from 
the charged Kerr metric, that the twist leads to intrinsic angular mo- 
mentum. The present equations of motion should, we believe, general- 
ize to motion for spinning particles. 

• We hope to generalize the second (or planned third) order calculations 
given here for the RTM equations to the general case of asymptotically 
fiat Einstein-Maxwell metrics. The condition of being algebraically 
special would be replaced by the tetrad condition of finding and using 
null geodesic congruences that are asymptotically shear-free[13, 15, 9]. 
To do this appears to be a quite formidable task. 
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8 Appendix A 

For completeness, though we have not needed them, the Maxwell field up to 
second order with / = 2 terms is given by 





q[2v 




Ife 




336\/2 
7 



with each expression being a function of r. 
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